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Abstract: The utilization of remotely sensed observations for light use efficiency (LUE) 
and tower-based gross primary production (GPP) estimates was studied in a USDA 
cornfield. Nadir hyperspectral reflectance measurements were acquired at canopy level 
during a collaborative field campaign conducted in four growing seasons. The 
Photochemical Reflectance Index (PRI) and solar induced chlorophyll fluorescence (SIF), 
were derived. SIF retrievals were accomplished in the two telluric atmospheric oxygen 
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absorption features centered at 688 nm (O 2 -B) and 760 nm (O 2 -A). The PRI and SIF 
were examined in conjunction with GPP and LUE determined by flux tower-based 
measurements. All of these fluxes, environmental variables, and the PRI and SIF exhibited 
diurnal as well as day-to-day dynamics across the four growing seasons. Consistent with 
previous studies, the PRI was shown to be related to LUE (r 2 = 0.54 with a logarithm fit), 
but the relationship varied each year. By combining the PRI and SIF in a linear 
regression model, stronger performances for GPP estimation were obtained. The strongest 
relationship (r 2 = 0.80, RMSE = 0.186 mg CCVnrVs) was achieved when using the PRI 
and SIF retrievals at 688 nm. Cross-validation approaches were utilized to demonstrate 
the robustness and consistency of the performance. This study highlights a GPP retrieval 
method based entirely on hyperspectral remote sensing observations. 

Keywords: gross primary production; light use efficiency; photochemical reflectance 
index; solar induced fluorescence; cornfield 


1. Introduction 

Carbon sequestration by terrestrial ecosystems is a key factor underpinning a comprehensive 
understanding of the carbon budget at a global scale. Terrestrial plants fix carbon through 
photosynthesis, known as gross primary production (GPP) at the ecosystem scale — the largest carbon 
exchange between the biosphere and the atmosphere [1]. Accurate measurements and estimates of GPP 
will allow us to achieve improved carbon monitoring and to quantitatively assess impacts from climate 
changes and human activities [2,3]. Remote sensing observations provide a unique opportunity to track 
photosynthetic activities at various spatial and temporal scales, and much effort has been put towards 
this goal over the past decades. Many remote-sensing-based indices and algorithms were developed 
and utilized to track leaf biochemical properties (e.g., chlorophyll, water content) and canopy 
biophysical properties (e.g., leaf area index, LAI; absorbed photosynthetically active radiation, APAR; 
fraction of APAR, fAPAR) have been applied to the estimation of GPP in various ecosystems, with 
variable success [2,4-14]. For instance, some studies reported good performances for GPP estimation 
(r 2 = 0.67 to 0.95) for agricultural sites [6,14] and grasslands [13] while less than satisfactory results 
(r 2 < 0.48) have been found for hardwood forests [13]. 

The light use efficiency (LUE) model [15,16] has been widely utilized and coupled with remote 
sensing data to estimate GPP at various spatial scales [17-22], including the Moderate Resolution 
Imaging Spectroradiometer (MODIS) Terra/Aqua satellite GPP Products [5,23,24]. The LUE concept 
measures the ability of vegetation to utilize atmospheric carbon dioxide (CO 2 ) and solar energy to 
produce biomass, and is usually defined as the ratio of GPP to APAR over a given period of time: 

GPP 

LUE = PAR x fAPAR (1) 

Whether based on biomass accumulation over foil growing seasons as originally formulated or as 
instantaneous flux measurements, LUE is highly variable among different species and over phenological 
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and environmental conditions. Instantaneous estimates of LUE are especially vulnerable to large 
uncertainties when used in GPP modeling [3,18,21,22,25,26]. For simplicity, the maximum expected 
LUE (LUE max ) has often been used to model GPP (e.g., MODIS GPP Product), but this approach 
requires considerable environmental information to down-scale from assumed optimal values to predict 
actual ones [23,24] and to describe GPP under non-optimal, unfavorable environmental conditions. 

LUE has been shown to be closely correlated to biochemical responses related to photoprotection of 
plant tissues, when more solar radiation is absorbed than can be utilized for photosynthesis [27]. One 
such process is non-photochemical quenching (NPQ), regulated by the xanthophyll cycle, in which the 
pigment violaxanthin is reversibly de-epoxidized to zeaxanthin via antheraxanthin [18,27-33]. The 
xanthophyll cycle can be monitored with the Photochemical Reflectance Index (PRI), a normalized 
difference band ratio based on the narrow spectral bands centered at a physiologically active 531 nm 
response and a reference band insensitive to the xanthophyll signal (e.g., 570 nm), as proposed by 
Gamon and colleagues [32,34]. The capability of the PRI to track photosynthetic activities, including 
LUE, has been investigated across plant functional types at leaf, canopy, and landscape levels 
using various instruments [19,20,29,32,34-41]. Moreover, many recent studies have investigated 
various factors that influence the PRELUE relationship at canopy or ecosystem scales including 
viewing geometry, canopy structure, leaf area index (LAI), soil background, pigment content 
(e.g., carotenoids/chlorophyll ratio), and shadow fraction [17,18,21,22,42-51], 

Another important photoprotection process is chlorophyll fluorescence, where excessive energy is 
expelled in order to prevent harmful photo-oxidation [3,18,28,32], Decades of laboratory studies and 
recent field-based solar-induced fluorescence (SIF) studies have shown that fluorescence provides a 
direct indicator of plant photosynthetic function (e.g., carbon fixation), enabling early detection of 
environmentally induced stress [52-60]. However, retrieving the SIF signal from passive remote 
sensing spectral observations over vegetation canopies under ambient solar illumination has been 
challenging because it is a weak signal added to the reflected radiance. Nevertheless, recent studies 
have shown success for SIF retrievals at leaf, canopy, and even satellite levels. Previously, since the 
SIF signal is small and requires very high spectral resolution for retrieval, attempts were made to 
develop optical reflectance-based approaches for SIF [61-65], In contrast to these relative indices, 
radiance-based algorithms enable direct estimates of SIF in physical or arbitrary units [54,55,66-70]. 

The chlorophyll fluorescence emission spectrum occurs across the red and near-infrared spectral 
region (650-800 nm) with two peaks near 685 nm and 740 nm. Radiance-based SIF retrieval 
approaches utilize measurements within and around atmospheric Fraunhofer lines that overlay the 
emission window to decouple SIF from reflected radiance. This is possible because the SIF signal 
contributes more to reflected radiance from vegetation canopies within these narrow dark lines where 
irradiance is low because of strong atmospheric absorption [54,55,57], The two telluric oxygen bands, 
O 2 -A (centered at 760 nm, ~1 nm width) and O 2 -B (688 nm, ~4 nm width) bands, are the most 
exploited absorption features for SIF retrievals. Many studies have reported that SIF can be detected 
using this passive remote sensing based concept [26,54,69,71-75]. The successful SIF retrievals to 
date have utilized the far-red (O 2 -A) retrieval to improve estimates and/or reduce uncertainty of 
photosynthetic activities, LUE and GPP, at leaf [76] and canopy levels [26,60,73,74,77]. However, 
few studies have successfully retrieved red fluorescence from the O 2 -B feature, nor have they 
examined red fluorescence in conjunction with canopy GPP measurements. 


Remote Sens. 2013, 5 


6860 


The fluorescence emission and NPQ energy dissipation are important components in the carbon 
fixation machinery of plants. The two processes are considered to be photosynthetic stress indices 
under sub-optimal environmental conditions and direct probes of plants’ physiological and 
photosynthetic status [29,32,53-56,78-81], In recent years, advances in remote sensing technology, 
especially those providing higher spectral resolution and signal-to-noise ratio, have enabled retrievals 
from high spectral resolution optical observations that correlate to the two dissipative pathways. Since 
both NPQ and chlorophyll fluorescence provide information relevant to photosynthetic activities 
(e.g., LUE and GPP) at the canopy scale [54,82,83], further investigations are warranted. In particular, 
the relationships among photosynthesis, PRI, and SIF derived from remote sensing observations, need 
further investigation, including examination of the red SIF. Several recent studies investigated the use 
of PRI and SIF to model FUE and GPP with various success [26,73,74,78], but the mechanism linking 
PRI and SIF to each other and to GPP is not yet well-described. In the current study, we introduce a 
collaborative field campaign conducted in a cornfield during four years. Hyperspectral observations 
and eddy covariance measurements were acquired during each year of the campaign. The study aims to 
investigate a greater range, both diurnal and seasonal changes, than previous publications. The study 
was undertaken to examine the capability of PRI and both far-red and red SIF to model the diurnal 
variations of two important photosynthetic parameters, FUE and GPP for multiple growing seasons. 

2. Methods 

2.1. Study Site and Field Data Collection 

The study site is located within the Optimizing Production Inputs for Economic and Environmental 
Enhancement (OPE3) experimental field (39.030°N, 76.845°W) maintained by the USDA Beltsville 
Agricultural Research Center (BARC) in Maryland, USA. Field campaigns were conducted on corn 
( Zea mays F.); different com cultivars were planted during each of the four growing seasons in 2008, 
2010, 2011, and 2012 (Table 1). Measurements were made on a total of 30 sunny days spanning 
different growth stages, ranging from early vegetative stages when the corn crop was actively growing 
at ~1 m tall and had produced 7-9 leaves (e.g., V7, V9) to fully mature canopies in the reproductive 
stage at -2 m tall (e.g., VT, Rl), through advanced reproductive development in the senescent stage 
(e.g., R4, R6). On each field day, measurements were taken multiple times between 9 am to 6 pm local 
time with approximately one-hour intervals, to sample the diurnal course. Sampling was done every 
one meter along a marked 1 00 m north-south direction transect in the middle of the field to minimize 
disturbance to the field. One complete set of sampling the transect was usually done within 30 min. 
Average of the samples was calculated to represent the field for the specific sampling time [84]. 

Canopy level hyperspectral reflectance spectra (400-1,000 nm, ~1.5 nm Full Width Half 
Maximum; FWHM) were obtained using a USB4000 Miniature Fiber Optic Spectrometer (Ocean 
Optics Inc., Dunedin, FF, USA) with a bare fiber. All spectral observations were acquired at nadir, 
above the canopy at a height of approximately 1 m. This was accomplished by placing the fiber optics 
on a height-adjustable pole-mount, where a custom-made fixture was designed to position the 
instrument at a desired view zenith angle and relative azimuth angle [85]. Incident solar irradiance was 
determined using a Spectralon reference panel (Fabsphere, North Sutton, NH, USA). Incident solar 
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irradiance was measured in between approximately every five samples. An additional spectrometer 
was set up to measure the reference panel continuously to monitor detailed changes of irradiance 
during the field days. Supplemental measurements included LAI and fPAR. LAI was measured with 
the LI-COR LAI-2000 plant canopy analyzer (LI-COR, Lincoln, NE, USA) and IPAR was determined 
using the LI-COR LI-191 line quantum sensor (LI-COR, Lincoln, NE, USA). 

Table 1 . Information about crop and environmental conditions in the USDA experimental 
cornfield in the four growing seasons. Dates are in the format of day of year (DOY). Total 
precipitation and average temperature are for the time period in Figure 1. Gross primary 
production (GPP) is in the unit of g C02/m 2 /d. 


Year 

Planting Date 

Varieties 

Maximum LAI 

Maximum GPP; 
Date 

Total Precipitation 
(mm) 

Average Temperature 
(°C) 

2008 

180 

TA 560-00 

3.27 

73.98; 214 

256.54 

20.80 

2010 

136 

Pioneer 35F37 

2.48 

75.17; 199 

410.22 

23.22 

2011 

145 

Pioneer 35K09 

2.79 

60.81; 192 

354.08 

24.51 

2012 

138 

Dekalb 57-67 

3.42 

59.66; 208 

291.09 

22.97 


2.2. Spectral Data Processing 


The PRI is usually calculated using two green bands at 531 nm and 570 nm as a normalized 
index [32,34]: 


PRI = 


p531 — p570 
p531 + p570 


( 2 ) 


In this study, we used canopy reflectance acquired over corn canopies centered at 531.01 nm and 
570.08 nm at their native FWHM of 1.5 nm to calculate the PRI. The O 2 -A (760 nm) and O 2 -B 
(688 nm) absorption features were utilized to derive SIF signals in this study. The fluorescence 
retrievals were determined using a modified Fraunhofer line depth (FED) algorithm [66] on the 
spectral observations over corn canopies. The FED principle [86] utilizes incident solar irradiance and 
radiance measurements reflected off vegetation canopies. By comparing the measurements inside and 
outside (reference band) an absorption feature, the estimate of fluorescence is derived. The modified 
FED algorithm [66] utilizes measurements on both shoulders (instead of one, the left shoulder) 
of an absorption feature to construct the reference band for a more realistic description [54,70]. 
Detailed description and review of different algorithms to retrieve SIF can be found in previous 
publications [54,67,70], In this study, we utilized the bands centered at 685.57 nm, 687.11nm, and 
691.77 nm for red fluorescence retrievals and the bands at 757.86 nm, 760.86 nm, and 772.67 nm for 
far-red fluorescence retrievals. These bands are at the native sampling resolution of the instrument of 
approximately 0.2 nm. We use SIF (red) to denote retrievals derived in the O 2 -B feature centered at 
688 nm and SIF (far-red) for retrievals within the O 2 -A feature centered at 760 nm. Since values 
for SIF retrievals are also affected by the intensity of irradiance, we also derived SIF yield, a 
dimensionless quantity that is independent of the light level. SIF yield was calculated following 
methods developed in a previous study [57] as the ratio of SIF retrievals and the measurement outside 
the absorption feature (reference band) from irradiance, denoted as SIF (red) yield in the O 2 -B band 
region and SIF (far-red) yield for the O 2 -A feature. Normalized Difference Vegetation Index (NDVI) 


Remote Sens. 2013, 5 


6862 


was calculated using red (620-670 nm) and near-infrared (841-876 nm) bands similar to the spectral 
characteristics of MODIS land surface reflectance. 

2.3. Flux Data, LUE and GPP Modeling 

The cornfield is equipped with a 10 m high eddy-covariance-instrumented tower that measures 
fluxes from a field approximately 16 hectares in area. Our sampling was approximately 65 m away 
from the flux tower and within the footprint. Flux data are reported at 30 min intervals throughout 
all 24-h periods of the full growing season. Net Ecosystem Production (NEP) is determined from 
the measured net CO 2 exchange using the eddy covariance method, supported by the tower’s 
micro-meteorology measurements (e.g., PAR, temperature, wind speed and direction, precipitation). 
Data were processed following the FluxNet/AmeriFlux guidelines [87], and the NEP was partitioned 
into respiration and GPP [88], compatible with AmeriFlux standards. Air temperature was measured at 
four meters above the ground on the flux tower and reported at half-hourly intervals. Precipitation 
was measured with a tipping rain gauge mounted one meter off the ground at a USDA/BARC 
meteorological station. 

Corresponding flux data were extracted when spectral observations were made on each field day. A 
combined flux and optical dataset was used to examine the potential of modeling LUE and GPP as the 
product of PRI, SIF (red and far-red), or both optical variables. This goal was accomplished through 
linear regression analyses with the R language and software by JMP (SAS Institute Inc., Cary, NC, 
USA). The general model formulations are listed in Table 2. 

Table 2. General formulation of the light use efficiency (LUE) and GPP flux parameters 
examined in this study, developed using the Photochemical Reflectance Index (PRI) and 
solar induced chlorophyll fluorescence (SIF) optical variables acquired in the USDA 
Beltsville Agricultural Research Center (BARC) cornfield. Note: a, b, c, d, represent 
coefficients of the linear regression models from statistical analyses. 


Output Variable 

Predictor Variable 

Formula 

LUE 

PRI 

LUE = a + b x PRI 


SIF 

LUE = a + b x SIF 


PRI, SIF 

LUE = a + b x PRI + c x SIF + d x PRI x SIF 

GPP 

PRI 

GPP = a + b x PRI 


SIF 

GPP = a + b x SIF 


PRI, SIF 

GPP = a + b x PRI + c x SIF + d x PRI x SIF 


2.4. Cross-Validation 

To evaluate the performance of the experimental models, cross-validation procedures were utilized. 
First we used the k-fold cross-validation approach. This method partitions the data into k subsets, then 
uses (k - 1) subsets as the training set to fit the model while the validation is conducted using the 
omitted subset. In this study, the k value was set to 10 and the dataset was randomly divided into 
10 subsets of equal sample size. The process was repeated 10 times to assess the coefficient r 2 cv and 
RMSEcv for performance evaluation. Secondly, we partitioned the dataset in a more restrictive manner, 
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utilizing the 2008 data (roughly half of the entire dataset) as the training data versus the 2010 through 
2012 data as the validation data. This strategy was invoked to assess the consistency of the model in 
the temporal domain. 

3. Results 

3.1. Diurnal and Seasonal Courses of GPP, PR1, and SIF 

The different cultivars planted each year produced variable leaf area indices (LAI) and maximum 
mid-season GPP, having been influenced by planting date and environmental conditions (Table 1). 
Daily GPP derived from tower-based eddy covariance measurements in the cornfield over four 
growing seasons is shown in Figure la. The grey dashed lines indicate the days when field campaigns 
were conducted. Corresponding daily precipitation, daily PAR, and daily average and high temperature 
observations are reported in Figure lb-d, respectively. Clear seasonal dynamics were observed for 
GPP, which increased during the green-up vegetative stage, reached the highest productivity during 
mid-season, and steadily declined after that (Figure la). Within this general seasonal trend, an increase 
in GPP was usually observed after a noticeable precipitation event (Figure la,b). The highest mid-season 
GPP value (~75 g CCVnT/d) was observed in early August of 2008 and in mid-July of 2010 
(Figure la; Table 1), and comparable GPP values and seasonal patterns were observed in those two 
years. Lower maximum mid-season GPP values were observed in 2011 and 2012, mostly due to 
unfavorable environmental conditions, such as unseasonably high temperature and/or lower precipitation 
(Figure lb,d). 

Figure 1. Seasonal variations of (a) daily GPP (g CCL/nTVd); (b) daily precipitation (mm); 

(c) daily PAR (mol/m 2 /d); (d) daily average and high air temperature in USDA BARC 
OPE3 cornfield for four growing seasons. X-axis values indicate day of year (DOY) in 
2008, 2010, 2011, 2012. The vertical dashed grey lines in panel (a) indicate field days. 
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Figure 1. Cont. 
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Figure 2. Diurnal patterns of (a) GPP, (b) LUE, (c) SIF (red), and (d) PRI for selected 
field days in four growing season. 
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Our dataset consists of a total of 30 days of carbon fluxes and spectral observations. The diurnal 
patterns of GPP, LUE, SIF (red), and PRI are shown in Figure 2a-d, respectively, for a sample dataset 
of eight days. In general, all four variables showed different maximum (or minimum) values and 
clearly expressed diurnal dynamics (Figure 2). Lower GPP values occurred in early morning and late 
afternoon, and higher values occurred around mid-day, in concert with illumination and air 
temperature changes (Figure 2a). In contrast, LUE exhibited higher values in early morning and late 
afternoon, with lower values during mid-day, indicating lower efficiency due to mid-day thermal or 
moisture stress (Figure 2b). SIF (red) and GPP exhibited similar diurnal courses (Figure 2c). Diurnal 


Remote Sens. 2013, 5 


6865 


changes of PRI and LUE showed similar patterns, with higher values in early morning and the late 
afternoon, but lower values during mid-day (Figure 2d). 

The seasonal changes of NDVI, PRI and SIF are shown as daily average values, along with daily 
GPP values, in Figure 3. NDVI expressed little seasonal change during the 2008 field season but had 
more variation in the other years. However, the PRI and SIF values exhibited pronounced variations in 
magnitude and seasonal dynamics across the four growing seasons. Towards the end of the growing 
season in 2008, 2011, 2012, the PRI values declined rapidly to notably lower values during the 
senescent stage (Figure 3). In general, the PRI and all of the SIF variables had seasonal trends similar 
to GPP. 


Figure 3. Cornfield seasonal dynamics of the daily average over four growing seasons for: 
(a) GPP, (b) PRI and Normalized Difference Vegetation Index (NDVI), (c) SIF yield 
(xl00%), and (d) SIF (W/m 2 /sr/nm). 
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3.2. LUE and GPP Modeling 

We investigated using PRI and SIF to model two important photo synthetic parameters, LUE and 
GPP. The statistics of modeling LUE using PRI and SIF are summarized in Table 3. In general, SIF by 
itself did not generate strong relationships with LUE when examined over time. By itself, however, the 
PRI exhibited the best correlation with LUE (r 2 = 0.45, RMSE = 0.000324 mg CCVpmol PAR). 
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Nevertheless, as shown in Figure 4, the relationship between PRI and LUE was not consistent among 
the four years (grey lines). Linear fit with r 2 of 0.56, 0.69, 0.50 was found for 2008, 2011, 2012, 
respectively with different slopes and offsets (Figure 4). By applying a logarithm fit, a 9% improvement 
was observed for the PRI model over the four-years (r 2 = 0.54, RMSE = 0.000322 mg CCVpmol PAR, 
Figure 4, black line). By combining the PRI with any SIF variable, substantial improvements in 
relationships to LUE and GPP were achieved. For LUE, a 16% improvement over the linear PRI 
model, and a 7% further improvement over the log PRI model, was obtained for the model that 
combined the PRI with the SIF (red) yield (r 2 = 0.61, RMSE = 0.000275 mg CCVpmol PAR; Table 3). 
The statistics of modeling GPP using PRI and SIF are summarized in Table 4. For GPP, >24% 
improvements over the linear PRI model (r 2 = 0.54, RMSE = 0.2770 mg C(Vm 2 /s) were obtained for 
the models that combined the PRI with either the SIF (far-red) (r 2 = 0.78, RMSE = 0.1894 mg C(Vm 2 /s) 
or with the SIF (red). The SIF (red) retrieval provided the strongest performance overall (r 2 = 0.80, 
RSME = 0.1994 mg CCVm 2 /s; Figure 5; Table 4). 


Table 3. Summary of statistics (r 2 and RMSE) of LUE models examined in this study. 
Note: the highest performing model in each group is shown in bold type. 


Output Variable 

Predictor Variable 

r 2 

RMSE (mg CCVpmol PAR) 

LUE 

PRI 

0.45 (0.54 Logarithm Fit) 

0.000324 (0.000322 Logarithm Fit) 


SIF (red) 

0.12 

0.000409 


SIF (far-red) 

0.01 

0.000435 


SIF (red) yield 

0.29 

0.000368 


SIF (far-red) yield 

0.06 

0.000424 


PRI, SIF (red) 

0.55 

0.000297 


PRI, SIF (far-red) 

0.48 

0.000317 


PRI, SIF (red) yield 

0.61 

0.000275 


PRI, SIF (far-red) yield 

0.53 

0.000301 


Figure 4. Relationship between canopy PRI and flux tower-based LUE for 2008 (♦), 2010 
(■), 2011 (A), 2012(X) in the USDA BARC cornfield. 
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Table 4 . Summary of statistics (r 2 and RMSE) of GPP models examined in this study. 
Note: the highest performing model in each group is shown in bold type. 


Output Variable 

Predictor Variable 

r 2 

RMSE (mg C0 2 /m 2 /s) 

GPP 

PRI 

0.54 

0.2770 


SIF (red) 

0.31 

0.3598 


SIF (far-red) 

0.28 

0.3891 


SIF (red) yield 

0.21 

0.3877 


SIF (far-red) yield 

0.20 

0.4107 


PRI, SIF (red) 

0.80 

0.1894 


PRI, SIF (far-red) 

0.78 

0.1994 


PRI, SIF (red) yield 

0.67 

0.2055 


PRI, SIF (far-red) yield 

0.66 

0.2099 


Figure 5. Performance of the GPP model (r 2 = 0.80, RMSE = 0.186 mg C02/m 2 /s) using 
both the PRI and the SIF (red) observations from four growing seasons (2008, 2010-2012) 
in the USDA BARC cornfield. 



Observed GPP (mg CO,/m 2 /s) 


3.3. Cross-Validation 

The performance of the GPP models, using cross-validation procedures, was examined for the model 
that combined the PRI with the SIF (red), since it produced the best empirical performance. The k-fold 
cross-validation approach to the model (Figure 6) consistently provided satisfactory statistics (r 2 cv = 0.79, 
RMSEcv = 0.188 mg CCVnr/s). When repeating this analysis for the same model using only 2008 data, 
evaluated by comparison with the remaining dataset, showed comparable results for the training data 
(Figure 7a) and validation data (Figure 7b). These results were also similar to the original model 
presented previously (Figure 5), demonstrating the consistency and robustness of the model across 
multiple growing seasons, including different corn cultivars and variable environmental conditions. 


Estimated GPP (mg C0 2 /m 2 /s) 
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Figure 6. Results of the k-fold cross-validation (. K value set to 10) performed on the 
GPP model presented in Figure 5. The cross-validated r 2 (r 2 cv = 0.79) and RMSE 
(RMSE CV = 0.188 mg CCVnr/s) confirms the robustness of the GPP model using PRI 
and SIF (red). 
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Figure 7. Cross-validation of the GPP model presented in Figure 5 using a systematic 
partitioning approach: (a) data acquired in 2008 were used as training data to develop the 
model; and (b) data from 2010 to 2012 were used as validation data. The results confirm 
the consistency of the GPP model applied over four growing seasons. 
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4. Discussion 

In this study, we examined two different types of optical remote sensing variables (PRI, SIF) to 
track the physiological status of plants for monitoring carbon assimilation in a cornfield. We made 
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observations over different environmental and crop conditions in four growing seasons, which 
included different corn cultivars in each year. 

We documented inter-annual and diurnal variations in the optical indices as well as GPP values 
(Figure 1), which were influenced by weather conditions. Precipitation in this rain-fed field increased 
GPP values. But unseasonably high temperatures, observed in the middle of the 2011 season and in 
2012 during the vegetative stage, created strongly expressed canopy stress which reduced productivity 
and GPP. Diurnal GPP cycles showed higher mid-day production when the vegetation was exposed to 
abundant solar radiation. Flowever, both the SIF and PRI responses indicated that plants experienced 
physiological stress during mid-day, also induced by high solar radiation and temperatures. Fligher SIF 
and lower PRI values expressed higher stress conditions during mid-day, while lower SIF/higher PRI 
values described stress abatement during early morning and late afternoon (Figure 2c, d). This diurnal 
trend was captured in the tower-based LUE, which exhibited the opposite diurnal dynamics to GPP 
(Figure 2a, b). A similar pattern to diurnal trends was seen in the seasonal trends of GPP, PRI, and SIF 
(Figure 3). Certain growth variables (e.g., LAI, phenological stage) affected the observed values for 
GPP, PRI and SIF. In 2008, NDVI value remained fairly consistent across the field season. On the 
contrary, both PRI and SIF displayed seasonal changes in concert with GPP that year, similar to the 
results for a sorghum field [75]. This shows that the PRI and SIF provide information that is more 
related to photosynthetic function than might be inferred from the NDVI. 

Several previous studies investigated similar topics and evaluated the performance of PRI and SIF 
for carbon monitoring. Rossini et al. [73] examined SIF (far-red) retrievals at 760 nm for fAPAR 
estimates, combined with either SIF yield or PRI, to model LUE in a rice field during two growing 
seasons; this method significantly improved estimates and reduced uncertainty in mid-day GPP 
(r 2 = 0.91, RMSE = 3.40 pmol CCh/nr/s). Rossini et al. [74] examined various formulations that 
coupled a suite of spectral vegetation indices, including the PRI, to estimate mid-day and daily GPP in 
a subalpine grassland. Damm et al. [26] utilized SIF retrievals at 760 nm and PRI to improve LUE 
estimates for describing the diurnal GPP course in a cornfield. Zarco-Tejada et al. [79] reported that the 
PRI and SIF retrievals at the 760 nm showed a similar seasonal trend to that of GPP measured at the same 
time of spectral observations (r 2 between 0.75 and 0.84; p < 0.01). The study reported a GPP model 
(r 2 = 0.67, RMSE = 5.77 pmol CCL/nr/s) using LUE estimates based on the SIF (far-red) yield at 760 nm, 
while the PRI-based LUE exhibited weaker performance (r 2 = 0.26, RMSE = 9.25 pmol C02/m 2 /s) for 
GPP estimates. These field studies demonstrate that variable results across vegetation types and 
conditions can be obtained in using SIF and/or PRI to estimate LUE or GPP. 

Our study addressed a similar topic: to investigate the potential use of PRI and SIF to estimate GPP, 
and to reduce uncertainties in the GPP and LUE estimates (Figure 5, Table 4). The models investigated 
and presented in our current study were composed solely of remotely-sensed observations. Unlike 
many previous studies, no micrometeorological variables (e.g., PAR or air temperature) were needed 
to achieve acceptable performance. The format of the models was determined after reviewing previous 
studies [26,73,74]. The models we utilized are similar to the Type 1 model in Rossini et al. [74], or the 
E3 model examined in Rossini et al. [73]. This confirms the potential for this type of GPP modeling 
approach, which has the advantage of simplifying the input data type requirements. Although most 
previous studies have focused on retrievals using the O 2 -A band to obtain SIF (far-red), we included 
SIF retrievals at the O 2 -B band feature to obtain SIF (red), which provided the strongest relationship to 
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GPP, possibly due to its important role in photosystem II (PSII) processes [53,89]. Furthermore, we 
undertook the challenge to model the diurnal variation of GPP during four growing seasons, which 
presented a complex task since both plant physiological conditions and GPP were changing 
significantly within any day, as well as across the four growing seasons due to varying climate 
conditions (Figure 1). 

The results (Tables 3 and 4) showed the challenge of modeling diurnal variations of LUE and GPP 
at this site for multiple growing seasons when using only PRI or SIF. Better correlations between PRI, 
SIF, LUE and GPP were found when using daily average values to study seasonal variations. For 
instance, for the 2008 growing season, when using daily average values, better performance was 
observed for the PRELUE model (r 2 = 0.76), the SIF (red) yield:LUE model (r 2 = 0.71), and the 
SIF (red):GPP model (r 2 = 0.61). However, the performance changed when we made the attempt to 
model diurnal variations across multiple growing seasons. In general, the GPP models examined 
in this study performed better than the LUE models. Multiple studies have demonstrated the 
potential of using the PRI to track LUE and photosynthetic activities for various vegetation types 
(e.g., [17,18,21,22,29,50,51]). The PRELUE model we developed with our dataset performed 
adequately (r 2 = 0.54) and overall results were comparable to previously published results [50], 
Nevertheless, one can easily observe year-to-year variations among the four growing seasons, since 
many factors affect the PRI values as well as the PRELUE relationship, e.g., canopy structure, LAI, 
soil background, and pigment content [42,45,48,50,51,85], No significant correlations were observed 
between SIF and either LUE or GPP for this dataset, but adding SIF to a simple PRI-based model 
improved its capability to track the variation in LUE by as much as 16% (Table 3). These findings 
suggest a future approach to significantly improve the performance of current LUE models, by 
combining these two optical signals. 

NPQ and chlorophyll fluorescence emission are considered as direct indicators of 
photosynthetic activities and physiological stress when plants are under sub-optimal environmental 
conditions [28,53,59,80,90-93] such as low water availability [81,93-95], unfavorable 
temperatures [96,97], nutrients [98-101], and salinity [102,103], Consequently, many studies have 
utilized the spectral indices, PRI and SIF retrievals, as indicators of stress and down-regulation of LUE 
and GPP [29,50,51,65,82,92,104-107], Our results lend support to those reported from previous 
studies that demonstrated that PRI and SIF retrievals were correlated to two important protective 
mechanisms (NPQ and chlorophyll fluorescence emissions), discarding excess absorbed solar radiation 
that could not be used in photosynthesis [3,29,54]. These might explain the empirical GPP model we 
presented in this study using both PRI and SIF showed good agreement with tower-based fluxes with 
low uncertainties. Furthermore, we achieved significant improvement by using both PRI and SIF 
together for GPP estimates, rather than using either PRI or SIF alone to attain a more comprehensive 
description of the entire carbon fixation mechanism, and hence, to deliver more accurate 
GPP estimates. 

Cross-validation approaches are designed to evaluate the robustness and consistency of statistical 
models. One crucial factor about cross-validation is the strategy to divide the entire dataset into 
training and validation data, using a random sampling and a systematic partitioning to test our results. 
We demonstrate that the GPP model based on 2008 data was able to predict GPP using PRI and 
SIF (red) in the three other years — 2010, 2011, and 2012. This highlights the robustness of the 
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integration of PRI and SIF (red) to capture important information about physiological status of this 
cornfield, thus providing accurate estimates of carbon assimilation with small uncertainties. The k-fold 
cross-validation approach was also applied to the model with PRI and SIF (far-red), with only a 
slightly lower performance (r 2 cv = 0.77, RMSE CV = 0.205 mg CCVnr/s). It can be difficult to retrieve 
SIF (red) from space [108] since it is a narrower and weaker feature than SIF (far-red), and this 
capability has yet to be attained. Even though our results suggested better performance of SIF (red) to 
model photosynthetic processes in ground level studies, we also demonstrated the potential for using 
SIF (far-red) for carbon monitoring using satellite measurements at global level. 

5. Conclusions 

In this study, we examined the capability of the PRI and SIF to track plant photosynthetic function 
in order to estimate diurnal variations of GPP in a cornfield over multiple growing seasons. Our results 
showed that the PRI and the SIF successfully captured the diurnal and seasonal dynamics in plant 
physiological status, in response to environmental and phenological conditions. Similar diurnal and 
seasonal patterns were observed for PRI and tower-based LUE as well as SIF and GPP. Although the 
PRI delivered the best performance (r 2 = 0.54; RMSE = 0.000322 mg C02/pmol PAR) when using a 
single optical variable for modeling canopy LUE, it nevertheless, was inconsistent when applied 
among multiple years. The most significant result was obtained by integrating both the PRI and the 
SIF, especially the SIF (red), in an empirical model to monitor diurnal variations in GPP for four 
growing seasons (r 2 = 0.80; RMSE = 0.1894 mg CCh/nr/s), which showed a minimum of 26% 
performance improvement than using either PRI or SIF alone. The robustness and consistency of 
which was demonstrated by cross-validation approaches (r 2 cv = 0.79; RMSE CV = 0.188 mg CCVnr/s). 
These results demonstrate that the PRI and SIF, utilized together, provide a powerful tool for carbon 
monitoring. This concept must be pursued over various vegetation types, given the success we report 
here for a cornfield, but we recommend more emphasis on SIF (red) retrievals in the future. 
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